Chiral phase transition in the presence of spinodal decomposition 
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The thermodynamics of a first order chiral phase transition is considered in the presence of 
spinodal phase separation within the Nambu-Jona-Lasinio (NJL) modeL The properties of the 
basic thermodynamic observables in the coexistence phase are discussed for zero and non-zero quark 
masses. We focus on observables that probe the chiral phase transition. In particular, the behavior 
of the specific heat and entropy as well as charge fluctuations are calculated and analyzed. We 
show that the specific heat and charge susceptibilities diverge at the isothermal spinodal lines. 
We determine the scaling behavior and compute the critical exponent 7 of the net quark number 
susceptibility at the isothermal spinodal lines within the NJL model and the Ginsburg-Landau 
theory. We show that in the chiral limit the critical exponent 7 = f/2 at the tricritical point as 
well as along the isothermal spinodal lines. On the other hand, for finite quark masses the critical 
exponent at the spinodal lines, 7 — 1/2, differs from that at the critical end point, 7 — 2/3, 
indicating a change in the universality class. These results are independent of the particular choice 
of the chiral Lagrangian and should be common for all mean field approaches. 

PACS numbers: 25.75.Nq, 24.60.Lz 



INTRODUCTION 



A first order phase transition is intimately linked with 
the existence of a convex anomaly in thermodynamic 
pressure [l| , which can be probed only in non-equilibrium 
states. There is an interval of energy and/or baryon 
number density where the derivative of pressure with 
respect to the volume of the system is positive. This 
anomalous behavior defines a region of instability in the 
(T, nq)-planc with Uq being the net quark number den- 
sity. The region is bounded by the spinodal lines, where 
the pressure derivative with respect to volume vanishes. 
The derivative taken at constant temperature or entropy 
defines the isothermal and isentropic spinodal lines, re- 
spectively. Spinodal decomposition is thought to play a 
dominant role in the dynamics of low energy nuclear col- 
lisions in the regime of the first order nuclear liquid-gas 
transition [H, Furthermore, the consequences of spin- 
odal decomposition in the connection with the chiral and 
deconfinement phase transitions in heavy ion collisions 
have been discussed in Refs. P, S, H, I3j SUl ■ It has been 
argued that in the region of phase coexistence, a spin- 
odal phase separation leads to a significant enhancement 
of baryon 3 a-iid strangeness fluctuations [1]. Studies 
of non-equilibrium dynamics in the linear sigma model 
have shown that significant density inhomogeneities may 
be dynamically induced in a first-order chiral phase tran- 
sition 0. In a recent paper it was shown that the spin- 
odal instability is connected with divergent fluctuations 
of the net quark number, (NqNq) - (^g)(A^q) = VTxq, 
along the isothermal spinodal lines [8|. Consequently, 
the quark number susceptibility Xq diverges not only at 
the critical end point (CEP) in the QCD phase diagram 
but also at the (non-cciuilibrium) first order phase tran- 



sition. The close relation of the singularities at the CEP 
with those at the spinodal lines is a consequence of the 
thermodynamic relations connecting the compressibility 
of the system with the quark number susceptibility. In 
this paper we explore the critical behavior and properties 
of the fluctuations of conserved charges and the specific 
heat in the presence of spinodal phase decomposition us- 
ing the Nambu-Jona-Lasinio (NJL) model. In particu- 
lar, we study the influence of an explicit chiral symmetry 
breaking through the finite value of the quark mass on 
the thermodynamics across the ofF-equilibrium first or- 
der transition. The critical exponents of Xq the CEP 
and along the spinodal lines are computed in the NJL 
model and in the Ginzburg-Landau effective theory. In 
the chiral limit, the critical exponents at the CEP and at 
the spinodal lines are identical, while for finite vales of 
the quark mass, the exponents differ. Our calculations 
are done within the mean-field approximation to the NJL 
model. However, the results are generic and should ap- 
ply to all QCD-like chiral models. We stress at this point 
that for non-universal quantities, like the position of the 
CEP, one can expect only qualitative agreement between 
different models. On the other hand, for universal quan- 
tities, like the critical exponents, all models in the same 
universality class yield the same result 

The paper is organized as follows: In Section [2] we 
explore the general features of the spinodal instabilities 
within the NJL model. In Section [3] we discuss the crit- 
ical behavior of the quark number susceptibility in the 



The critical exponents obtained in the mean- field approximation 
arc in general modified by fluctuations. 
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mean field approximation. The critical exponents of Xq 
are computed in Ginzburg-Landau theory in Section 31 
Concluding remarks are presented in Section [5l 



2. EFFECTIVE CHIRAL MODEL AND 
SPINODAL INSTABILITIES 

We study the thermodynamics of the chiral phase tran- 
sition in the presence of spinodal instabilities within the 
Nambu-Jona-Lasinio (NJL) model [§]. The model de- 
scribes the effective interactions of quarks preserving the 
chiral symmetry of the massless QCD Lagrangian. We 
include quark interactions only in the scalar and isoscalar 
sectors. In this case the NJL Lagrangian for two quark 
flavors, three colors and for finite quark chemical poten- 
tials has the form [13, [HI, [H, d : 

(2.1) 

where m — diag(mti,TOd) is the current quark mass, 
fiq ~ diag(/i„,^d) is the quark chemical potential and 
T are Pauli matrices. The strength of the interactions 
among the constituent quarks is controlled by the cou- 
pling GsA^ = 2.44 with the three momentum cut-off 
A = 587.9 MeV, introduced to regulate the ultraviolet 
divergences [l3|- The parameters of the model are fixed 
so as to reproduce the pion mass and decay constant in 
vacuum for m„ = rud = 5.6 MeV. 

In the mean field approximation the thermodynamics 
of the NJL model, for an isospin symmetri c sy stem, is 
obtained from the thermodynamic potential [lO|: 
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where M ~ m — 2Gs{il'ip) is the dynamical quark mass, 
Eip) 



y/]p+M^ is its energy and n^^\p,T,^) 



1 + exp[(£^(p) =F /i)/r] j is the particle/antiparticle 

distribution function. The quark chemical potential /i is 
expressed as an average ^ = /i„ = fXd- The dynamical 
quark mass M in Eq. (|2.2[) is obtained self-consistently 
from the stationarity condition dQ,/dM = 0, which im- 
plies: 
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The thermodynamic pressure of the system is given 



by P = — ri/V^, and the net quark number density 
given by 
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Thus, in the mean-field approximation the net quark den- 
sity Uq has the same structure as in a non-interacting 
gas of massive fermions, albeit with a (T, /x)-dependent 
effective mass Al (T, /x) . The dynamical quark mass M 
is an order parameter of the chiral phase transition. In 
the left panel of Fig. [1] we show M computed at fixed 
temperature T = 30 MeV as a function of the quark 
chemical potential, in the chiral limit, i.e., for vanish- 
ing current quark masses. The behavior of M is typi- 
cal for an order parameter in a first order phase transi- 
tion. There is a range of values of the chemical potential 
where there are several solutions to the gap equation. 
There are stable/metastable solutions, corresponding to 
absolute/local minima of the thermodynamic potential 
and an unstable solution, corresponding to a maximum. 
The equilibrium transition from the chirally broken to the 
symmetric phase is obtained when the two minima are 
degenerate. This corresponds to the so called Maxwell 
construction. The location of the equilibrium transition 
is indicated in Fig.[l]by the dashed line. A non-zero cur- 
rent quark mass in the NJL Lagrangian explicitly breaks 
the chiral symmetry. Consequently, in this case, the dy- 
namical mass M is not a true order parameter and it 
never vanishes. The dependence of M on the quark chem- 
ical potential is shown in the right panel of Fig. [TJ The 
characteristic dependence of M on the chemical potential 
at a first order transition is, as will be discussed below, 
refiected in various thermodynamic quantities. 

In Fig. [2] we show the net quark number density at 
fixed temperature as a function of fx for vanishing as well 
as for a finite current quark mass. As expected, the non- 
monotonic dependence of the mass M on the chemical 
potential results in corresponding structure in Uq. How- 
ever, there is an notable difference in the properties of Uq 
in the chiral limit and for a finite quark masses. In the 
latter case, Uq is differentiable at all values of /i, while 
in the chiral limit Uq exhibits a cusp at the point where 
the dynamical quark mass vanishes. This has important 
consequences for observables obtained from Uq through 
differentiation with respect to thermal parameters, like 
e.g. the quark number or mixed susceptibilities. The 
thermodynamic pressure as a function of the volume per 
net quark, (l/uq), is shown in Fig.[3]at fixed temperature. 
Independently of the current quark mass, the qualitative 
structure of the pressure is that of a van-der-Waals' equa- 
tion of state. The unstable solution of the gap equation 
corresponds to the part of the equation of state where the 
volume derivative of the pressure is positive, i.e., to the 
mechanically unstable region. This region is bounded by 
the spinodal points located at the minimum and maxi- 
mum of the pressure, respectively. Outside of this region 
the system is mechanically stable "^'^ . The volume deriva- 



More precisely, the system is absolutely stable outside the coex- 
istence region, while within this region but outside the spinodal 
lines it is metastable, i.e., stable only against small amplitude 
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FIG. 1: The dynamical quark mass at fixed temperature (T = 30 MeV) as a function of the quark chemical potential in the 
chiral limit (left panel) and for a finite quark mass m = 5.6 MeV (right panel). The dashed line indicates the equilibrium first 
order phase transition while the dotted lines show the isothermal spinodal points. 
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FIG. 2: The net quark number density at T = 30 MeV computed in the chiral limit (left panel) and for a finite quark mass 
m = 5.6 MeV (right panel). The dashed and dotted lines have the same meaning as in Fig. [1] 



tive of the pressure depends on which quantity is kept 
constant, the temperature or the entropy. If the volume 
derivative of P exists, then the isothermal and isentropic 
spinodal lines are defined by zeroes of the corresponding 
derivatives 







dP 
dV 



0. 



(2.5) 



From Fig. [3l it is clear that for a finite quark mass m. the 
derivative of the thermodynamic pressure is well defined 
everywhere. However, in the chiral limit the pressure ex- 
hibits a cusp rather than a regular minimum, at the point 



fluctuations. 



where the pressure derivative changes sign on the high 
density side. Consequently, in this case equation (|2.5p 
can not be used to determine the position of the spinodal 
line. Instead, the corresponding branch of the isothermal 
or isentropic spinodal line is defined by a change in sign of 
the pressure derivative, or cquivalcntly by a minimum of 
the pressure. As wc discuss below, the corresponding sin- 
gularity of the thermodynamic potential is weaker than 
at a regular spinodal line. In Fig. [3] the phase diagram of 
the NJL model in the (T, nq) plane is shown, including 
the isentropic and isothermal spinodal lines. For finite 
m the NJL model yields a generic phase diagram with 
a first-order phase transition at low temperatures and 
high densities and a cross-over transition at high temper- 
atures. The critical end point separates the cross-over 
from the first-order transition. In the chiral limit, the 



4 




0.1625 



0.162 



0.1615 



0.161 



0.1605 



0.16 



0.1595 



1/nq [fm?] 




FIG. 3: In the left panel we show the pressure as a function of inverse quark number density for fixed temperature, T = 30 
MeV obtained in the chiral limit, while in the right panel we show the results for a finite quark mass, ?n — 5.6 MeV. 
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FIG. 4: The phase diagram of the NJL model in (T, n9)-plane. The left figure corresponds to the chiral limit. The dot indicates 
the position of the TCP and the full-line above the TCP shows the second order transition line. The full-lines below the TCP 
represent the phase boundaries of the first order chiral transition. Finally, the dashed (dotted) lines show the isothermal 
(isentropic) spinodal lines. The right figure was computed with a finite quark mass m = 5.6 MeV. Above the CEP, indicated 
by a dot, there is cross-over transition. 



cross-over transition becomes a second-order phase tran- 
sition. In an equilibrium first-order phase transition, the 
meta-stable and unstable regions are not probed. The 
states of maximum entropy, obtained by the Maxwell 
construction, correspond to the coexistence of the high- 
and low-density states at the corresponding phase bound- 
aries. On the other hand, in a non-equilibrium system, 
where the thermal parameters arc changing sufRcicntly 
fast, the system may traverse the meta-stable and enter 
the unstable region, bounded by the spinodal linos. From 
the thermodynamic relation 



it is clear that the isentropic spinodal lines are located 
inside the isothermal spinodal region and furthermore, 
that the two lines coincide at T = 0. As shown in Fig. 21 
the isothermal spinodal lines join at the CEP, whereas 
the boundary of the isentropic spinodal region remains 
well below the CEP. The latter property is probably an 
artefact of the mean-field approximation. When fluctua- 
tions arc included, the specific heat at constant volume 
Cv diverges at the CEP, although with a weaker singu- 
larity than Cp, while in the mean- field approximation 
Cv remains finite (see section 3). According to (|2.6p 
the isentropic and isothermal volume derivatives of the 
pressure arc equal for Cy oo. This implies that also 
dP/dV\s = at the CEP and that the isentropic spin- 
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odal region probably extends to the CEP. In order to 
explore this question in more detail and to determine 
the shape of the isentropic spinodal region, a system- 
atic study of fluctuations is required. At a first-order 
phase transition the instabilities are reflected in a con- 
vex structure of several thermodynamic functions. Such 
a structure is seen in the thermodynamic pressure and 
also appears in the entropy density as a function of quark 
density Uq 



-2N,Nf 
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However, the entropy per quark, s/uq^ exhibits a mono- 
tonic dependence on the density. In Fig. [5] we show s/nq 
computed at fixed temperature in the chiral limit and for 
a finite quark mass. In the chiral limit the entropy per 
quark shows a cusp structure at the spinodal line on the 
high-density side. 



SUSCEPTIBILITIES AND THE CRITICAL 
BEHAVIOR 



At the critical end point the fluctuations of the quark 
density diverge [H, [iB, [H, [13, E Ell , while they are fi- 
nite above and below the CEP if the first order transition 
takes place in equilibrium. However, in non-equilibrium 
the system may reach the spinodal instability. It is 
thus natural to explore the charge fluctuations in the 
mctastable and unstable regions of the first order phase 
transition. In statistical systems a measure of the fluc- 
tuations of a particular quantity is provided by the cor- 
responding susceptibility This are obtained as second 
derivative of the thermodynamic pressure. For the char- 
acterization of the spinodal region of the first order chiral 
phase transition, the susceptibilities related to the net 
quark number density Xq = X^^^J,l electric charge density 
Xq as well as the thermal xtt and mixed x^lT suscep- 
tibilities are of particular interest. These are obtained 
from the thermodynamic pressure (|2.2p as follows: 
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The susceptibilities involve derivatives of the dynamical 
quark mass with respect to fi or T Consequently, 
at the spinodal lines, where both dM/dfji and dM/dT 
diverge (see Fig. [Ij, one expects singularities in these 



'^■^ The above form of the mixed susceptibility is obtained from 
XiiT = ds/dfi. An alternative equation, involving dM/dT, is 
obtained from XpT = dnq/dT. The equivalence of the two is 
obvious when the explicit expressions for dAI/dfj, and dM/dT 
are inserted. 



susceptibilities. In Fig. [S] we show the quark number 
susceptibility as a function of Uq at a fixed temperature 
T = 30 McV. As expected, it is singular at the spinodal 
points. However, in the chiral limit the susceptibility 
exhibits only a weak singularity at the high-density point, 
corresponding to the cusp in the pressure in Fig. [3l At 
this point the susceptibility is discontinuous, changing 
from a finite negative to a finite positive value, while at 
the low density spinodal point the susceptibility diverges, 
corresponding to very large fluctuations of the net quark 
density. Moreover, for a finite current quark mass, the 



6 




0.5 I — ' ^ 1 1 1 1 1 • 1 0.5 I— ' 1 1 1 1 1 1 1 

0.2 0.4 0.6 0.8 1 1.2 1.4 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 

nq [fmr^] nq [fnr3] 



FIG. 5: The ratio of the entropy to the net quark number density, i.e., the entropy per quark, at fixed temperature as a 
fiinction of the net quark density in the chiral limit (left panel) and for a finite quark mass m = 5.6 MeV (right-panel). The 
dashed and dotted lines have the same meaning as in Fig. [1] 
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FIG. 6: The susceptibility at fixed temperature T = 30 MeV in the chiral limit (left panel) and for a finite quark mass 
m = 5.6 MeV (right panel). The vertical dashed and dotted lines have the same meaning as in Fig. [1] 



susceptibility diverges at both spinodal points. In both 
cases the quark number susceptibility is negative in the 
spinodal region, indicating a mechanical instability of the 
system. 



The critical properties of the net quark number sus- 
ceptibility seen at the isothermal spinodal lines are not 
there along the isentropic spinodal trajectories. There, 
the susceptibilities are finite and continues when the sys- 
tem passes through the metastable region. The above 
behavior of Xnn appears as a direct consequence of the 



thermodynamic relations 
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that connect the pressure derivatives with the suscepti- 
bilities Xxy = —d^n/dxdy. 

Along the isothermal spinodal lines the pressure 
derivative in Eq. (j3.4p vanishes. Thus, for non vanishing 
density Uq its fluctuations have to diverge. In addition, 
since the pressure derivative dP/dV\T changes its sign 
when crossing the spinodal trajectories, there is a corre- 
sponding sign change of the divergence in Xnn a-s seen in 
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Fig. [6l A similar behavior as for x^ifi is also expected for 
electric charge xq : strangeness xs or isovector xi fluctu- 
ations if their charge densities are finite. 

In heavy ion collisions the isospin asymmetry is usu- 
ally negligible at high energy and the strangeness den- 
sity vanishes due to strangeness neutrality in the ini- 
tial state. Thus, here not all chemical potentials are 
thermodynamically independent. The net quark num- 
ber Xq = —d'^il/dfig and the isovector xi = —d^il/d^'j 
susceptibilities are related with the electric charge fluc- 
tuations Xq as 



Xq 
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For isospin symmetric system the derivative of the pres- 
sure vanishes and all relevant susceptibilities are linearly 
dependent. Clearly, since Xg diverges on the isothermal 
spinodals and since n/ ~ is a good approximation in 
heavy ion collisions, thus the electric charge fluctuations 
XQ should show similar critical behavior to that seen in 
Xq- On the isentropic spinodals the quark number fluctu- 
ations arc finite. This is also valid for all mixed suscepti- 
bilities appearing in Eq. (j3.5p . The isentropic condition 
(|2.5p holds since the numerator in (|3.5p vanishes. The 
exception is the point at T = where the isothermal and 
isentropic spinodal conditions are equivalent since the en- 
tropy density and the Xp,T vanish at T ~ 0. The thermal 
Xtt and mixed XfiT susceptibilities along the isothermal 
spinodals show similar critical behavior as found for x^l^l 
in Fig. [HI When going beyond the mean field dynamics, 
the isothermal and the isentropic spinodal curves could 
also coincide at the CEP. Thus, from Eqs. and (|3?5l) 
it is clear that all susceptibilities should diverge at the 
CEP in this case. 

For equilibrium first order phase transition the prop- 
erties of charge fluctuations are different from those seen 
in Fig. [6l There, under the mean field approximation, 
the charge susceptibilities are positive and exhibit a dis- 
continuity at the transition point [lO, Hlj . This discon- 
tinuous structure is eventually converted to a cusp if the 
quantum correction are included [2l|, [l^]. However, the 
essential difference is that the susceptibilities in the equi- 
librium transition are finite whereas they are divergent if 
spinodal phase separation appears. 

The rate of change in entropy with respect to temper- 
ature at constant volume or pressure gives the specific 
heat expressed as 
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From Eqs. ([3T5|) . ((37f|) and ((3^ it is clear that the ratio 
of specific heats Cp/Cy satisfies a well known relation 
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In Fig. [7] we show the specific heat Cv at fixed temper- 
ature. One sees in the chiral limit a discontinuity at the 
right-handed spinodal branch which comes from a sim- 
ilar jump in Xxy At the left-handed spinodal branch 
Cv is continuous and finite. Each susceptibility Xxy di- 
verges at the spinodal point, however such singularities 
are canceled out in Eq. (|3.7p . The above is also valid 
for finite quark masses, thus Cy continuously changes at 
both spinodals. 

On the other hand, as seen in Fig. [3 the Cp diverges at 
the isothermal spinodals because there is no cancellation 
of the singularities among Xxy in Eq. (j3.8p . It has been 
argued that the negative specific heat could be a signal of 
the liquid-gas phase transition [2^ and its occurrence has 
recently been reported as the first experimental evidence 
for such an anomalous behavior in low-energy nuclear 
collisions [13, m. 

From Eqs. (|3.5p and (|3.8p one finds that Cp vanishes 
along the isentropic spinodal lines. However, this is not 
the case in the chiral limit. Here the pressure shows 
a cusp, thus its derivative is not defined. Consequently, 
dP/dV\s is non- vanishing at this isentropic spinodal line 
and Cp remains finite. 

At the CEP the fluctuations of charge densities that 
couple to the massless modes of the sigma field are di- 
vergent. Fig. [5] shows the x^^^ near the CEP and at fixed 
temperature corresponding to the first order transition. 
From this figure it is clear that when approaching the 
CEP from the side of the first order transition the re- 
gion of instabilities shrinks and at the CEP disappears 
completely. Consequently, the CEP singularity in 
appears from the matching of the two positive singular 
branches of Xmm- The strength of the singularity is con- 
trolled by the critical exponent 7 shown in Fig. [TOl In the 
chiral limit the exponent at the isothermal spinodal line 
is identical to that at the TCP. From numerical analysis 
we find that Xmm ~ ~ ^^c)~^ , with 7 = 1/2. How- 
ever, for finite quark masses we find 7 = 2/3 at CEP, in 
agreement with that expected in the mean-field approx- 
imation [21I, [2^ . From the numerical studies of Xmm at 
spinodals we found that 7=1/2 independently of the 
value of the quark mass. Thus, the singular properties 
at the TCP/CEP are the remnant of divergent behavior 
of the X/i/i at the first order chiral phase transition [1]. 
Clearly, the critical exponents are renormalized by quan- 
tum fluctuations [l^]. However, the smooth evolution 
of the singularity from the spinodal lines to the CEP is 
expected to be generic. 

In the next section we disuse how the critical exponents 
of quark susceptibilities along the spinodal lines can be 
analytically justified within the Ginzburg-Landau theory. 



8 




0.2 0.4 0.6 0.8 1 1.2 1.4 0.2 0.4 0.6 C.8 1 1.2 1.4 



nq [fm 3] nq [fm 3] 

FIG. 7: The specific heat Cv at constant volume and at fixed temperature T = 30 MeV in the chiral limit (left panel) and 
for a finite quark mass m = 5.6 MeV (right panel). The vertical dashed-lines indicate the location of the first order phase 
transition. The dotted-lines show the position of the isothermal spinodal points. 
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FIG. 8: The specific heat Cp at constant pressure and at fixed temperature T = 30 MeV in the chiral limit (left panel) and 
for a finite quark mass m = 5.6 MeV (right panel). The vertical dashed-lines indicate the location of the first order phase 
transition. The dotted/dashed-dotted lines show the position of the isothermal/isentropic spinodal points. 



4. GINZBURG-LANDAU EFFECTIVE THEORY 
AND CRITICAL EXPONENTS 

In the mean filed approximation and near the chi- 
ral phase transition the thermodynamic potential ob- 
tained in QCD-like chiral models can be described by the 
Ginzburg-Landau effective theory [2^. Consequently, 
the mean field critical exponents at the chiral phase 
boundary can be derived analytically. In earlier stud- 
ies within this approach the critical exponents of differ- 
ent susceptibilities were obtained assuming that the sys- 
tem appears in equilibrium [iqI . [2f| . In the following we 
extend the Ginzburg-Landau formalism to account for 
spinodal phase decomposition in the mcta-stable regime 
of the 1st order transition. 



According to the Ginzburg-Landau theory, close to the 
phase boundary, the thermodynamic potential may be 
expanded in power series of the order parameter M [26l |: 

n{T, M) = r^o - hM + -a2M^ + -a^W^ + -a^M^ , 

2 4 6 

(4.1) 

where hM is a term that explicitly breaks the symmetry 
AI ^ —M. In QCD this corresponds to the quark mass 
term. The gap equation 

n'{M) = -h + a2M + aiM^ + a^M^ = (4.2) 

yields the values of the order parameter that extrem- 
ize the thermodynamic potential. Consider first the case 
without explicit symmetry breaking, i.e. /i = 0. Then, 
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FIG. 9: The net quark number density fluctuations Xmm/^^ 
normalized to momentum cut-ofF as a function of tiie quark 
number density Uq. The full-line represents the fluctuations 
at the CEP at constant temperature Tcep — 80 MeV. The 
dashed-line shows the Xmm at flxed T — 70 MeV corresponding 
to the 1st order transition. 



At a spinodal line the gap equation is satisfied, i.e., 
Ai = 0. Thus, close to a spinodal line. 



Ai = ai (T - To) + bi{iJ - /io) 

A2 = a2(T - To) + 52(yU - yUo) , 



(4.7) 



while the coefficient A3 is non-zero, and can be taken 
to be constant over the small range of temperatures and 
chemical potentials considered. In the gap equation for 
a point close to the spinodal line /.i ~ ^0 (for definiteness 
we set T = To), 



Ai + A2m + Asm + Aiirr = , 



(4.8) 



we assume that m ~ ±|/z — /xo|", with < a < T For 
small l/i — /zo|, the leading terms of (|4.8p are then Ai + 
A^m? = 0, yielding a = 1/2. At the critical end point 
also vis = 0, so the leading terms in the gap equation 
are Ai + A^rm? ~ 0, and correspondingly a = 1/3. The 
singular part of the baryon number susceptibility is given 

by 



y 



b 



IM - Mol 



(4.9) 



at a second-order phase transition 02 = and 04 > 0, 
while at a first order transition 02 > 0, 04 < 0. The pa- 
rameter ag should be positive for stability. Clearly, the 
TCP can be identified by a2 = 04 = 0. If ft, 7^ and 
04 > 0, a change of sign of 02 corresponds to a cross over 
transition. We note that in this case it is not possible to 
uniquely define the position of the transition, since the 
order parameter M is always non-zero. The endpoint of 
the first order transition is then called a critical end point 
(CEP). In order to find the spinodal fines, we need also 
the second derivative of 



n"{M) = 02 + 804 A/^ + baeM^ . 

The zeroes of the second derivative, fi"(A/o) 
given by 



Af2 



3a4 



± 



1 



lOag lOag 



9al 



20a2ag . 



(4.3) 

0, are 

(4.4) 



Depending on the values of the parameters, the solutions 
correspond to inflexion points, or, if Mq also solves the 
gap equation, to a point on a spinodal line. We expand 
the thermodynamic potential around a zero of the second 
derivative 



A. 



A. 



n{m) = f^o + + — + —m^ + 
2 3 

where m = M — Mq. One finds 

Ai = -h + 02^/0 + 04^/0^ + aeM^ 



4 4 
m , 



(4.5) 



A2 = 02 
A3 = 80471/0 



8a4Mo^ + SaeAfo = 



9al - 20a2af{4.6) 



04 + IOgsMq = -2a4 ± 



9a| — 20a2aQ . 



with 7 = 1/2 at the spinodal lines and 7 = 2/8 at the 
CEP. Thus, the singularities of the two spinodal lines 
conspire to yield a somewhat different critical exponent 
at the CEP. In the chiral limit, a = 1/2 at the spinodal 
lines and the corresponding critical exponent is 7 = 1/2. 
Finally, at the TCP a = 1/4, but the critical exponent 
is unchanged, since here AI — and 



M- 



BM dM^ 



'1/2 



(4.10) 



Thus, in the chiral limit the critical exponents of the 
baryon number susceptibility at the spinodal lines and 
at the TCP are identical, while for finite quark masses 
the corresponding exponents differ. The critical and the 
universal properties at the phase boundary can be also 
studied considering the fluctuations of the order param- 
eter. They are characterized by the chiral susceptibility, 
Xmm defined as the second order derivative of thermody- 
namic potential with respect to the quark mass 



dmn? 



(4.11) 



Figure [TT] shows the change of the chiral susceptibility 
near the critical end point and along the spinodal lines 
in the chiral limit and for the finite value of the quark 
masses. In the chiral limit the behavior of Xmm as a func- 
tion oi t = |(// — /ic)/Mc| differs from that shown in Fig. 
[TPl for the net quark number density fluctuations Xnn- 



The assumption |a| > 1 docs not yield a consistent solution of 
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FIG. 10: The net quark number susceptibility in the vicinity of the TCP (left panel) and CEP (right panel) as a function 
of the reduced quark chemical potential t = {fi — fic)/f^c at fixed T. The filled circle (•) denotes the results calculated at 
T = Ttcp/cbp and the filled triangle (a) at T = 30 MeV < Ttcp/cep corresponding to the first order transition. 
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FIG. 11: The chiral susceptibility in the vicinity of the TCP and spinodals in the chiral limit (left panel) and at finite quark 
mass (right panel) as a function of the reduced quark chemical potential t — {fi — fic) / fJ-c at fixed T. The filled circle (•) denotes 
the results calculated at T = Ttcp/cep and the filled triangle (a) at T = 30 MeV < Ttcp/cep corresponding to the first order 
transition. 



First, the critical exponents at TCP and along spinodals 
do not coincide. Second, the critical exponent 7 at the 
TCP is now equal to unity instead of 7 = 1/2 whereas 
at spinodals it still remains the same as for X/^/i- On the 
other hand, at finite quark mass the critical exponents at 
the CEP and at the spinodals for Xmm are the same as 
found for Xmm - "^^^ above properties are to be expected 
in mean field dynamics because of the following relation 

X^^^ X7^' + M\mm, (4.12) 

that connects Xmm with Xmm and the dynamical quark 
mass M. The first term x^*^^ corresponds to a non- 
singular contribution to the quark number fluctuations. 
In the chiral limit, when approaching the TCP along any 



paths being not asymptotically tangential to the criti- 
cal line the order parameter ~ t^^'^ and at the 0(4) 
line ^ The corresponding exponents of Xmm 

arc 7 = 1, what explains through Eq. (|4.12p the critical 
exponents found in the Fig. [TTI-lcft. For the first-order 
transition in the chiral limit and for the finite quark mass 
the is always finite, consequently the critical expo- 
nents of Xfifi and Xmm are identical as seen in Fig. [TlJ 



5. SUMMARY AND CONCLUSIONS 

Considering deviations from an idealized equilibrium 
picture of the first order chiral phase transition by in- 
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eluding spinodal instabilities, we have explored the criti- 
cal behavior of the fluctuations of conserved charges and 
the specific heat in the NJL model in mean field approxi- 
mation. We have studied the properties these observables 
in the chiral limit and for finite vale of the quark chem- 
ical potential. We showed that the charge fluctuations 
diverge along the isothermal spinodal lines of the first 
order chiral phase transition independently on the value 
of quark masses. Accompanied by the singularities of 
the susceptibilities, the specific heat at constant pressure 
also diverges. This is in contrast to an equilibrium tran- 
sition where the susceptibilities are finite and only at the 
TCP/CEP diverge. This implies that large fluctuations 
of the density would be a signal not only for the exis- 
tence of the CEP or TCP but also of a non-equilibrium 
first order transition. The above is especially relevant 
to the current and future heavy-ion experiments when 
considering a possible signal for a first order chiral phase 
transition. 

We also showed that in the chiral limit the critical 
exponents of the net quark number susceptibility arc 
identical at the TCP and at the isothermal spinodals. 
In the mean field approximation the critical properties 
of the net quark number fluctuations were found to be 
governed by a common critical exponent 7 — 1/2. On 
the other hand at finite quark masses the exponents are 
different indicating a change of the universality class. At 
the spinodals the critical exponents for m ^ stays to 
be 7 = 1/2. However, at the CEP it changes from 7 = 
1/2 at m = to 7 = 2/3 for m 7^ 0. These values of 
the critical exponents were calculated numerically in the 
NJL model and also confirmed by the analytic studies of 
an effective Ginzburg-Landau theory. The above results 



indicates that the singularities of Xfifj, at the TCP and 
CEP are a remnant of divergent fluctuations along the 
two spinodal branches. 

The calculations were performed within the NJL model 
under mean field dynamics. However, the appearance 
of the singularities in the quark fluctuations and specific 
heat comes from the thermodynamic relations which con- 
nect the pressure derivatives with different charge suscep- 
tibilities. Thus, divergent properties of charge suscepti- 
bilities and the specific heat in the presence of spinodal 
instabilities are quite general, being independent on the 
specific formulation of the chiral Lagrangian. This is also 
the case for the fluctuation critical exponents at the spin- 
odals and their relation to the critical properties at TCP 
and CEP that should be valid for all QCD like chiral 
models under mean field approximation. It is a further 
challenge to verify the influence of the quantum and ther- 
mal fluctuations on the critical medium properties in the 
presence of spinodal phase separation. 
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